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The Tif = 2 residual mass in lattice HQET to order* 
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We compute the so called residual mass in Lattice Heavy Quark Effective Theory to order in the Uf — 2 
(unquenched) case. The control of this additive mass renormalization is crucial for the determination of the heavy 
quark mass from lattice simulations. We discuss the impact on an unquenched determination of the b-quark mass. 



1. Introduction 

The computation of quark masses is one of 
the most successful! application of Lattice QCD. 
Within this framework an interesting role is 
played by the determination of the 6-quark mass. 
Despite the fact that one can not accomodate 
the &-quark on a lattice, its mass is computed 
to a very good accuracy. As it is known, the key 
of the success is the use of some form of effec- 
tive theory. Many approaches have been used, 
whose results are consistent within errors in the 
quenched approximation [112131415] . A popular 
determination comes from Heavy Quark Effective 
Theory (HQET). It shares with other approaches 
the need of dealing with the mass counterterm 
we will be concerned with. Within HQET, the 
most direct relation between the mass of a physi- 
cal hadron {Mb) and the mass of the heavy quark 



Mb =mb + £ + 0{l/mb), 



(1) 



in which rrib is the HQET expansion mass pa- 
rameter and £ is the (linearly divergent) binding 
energy. However, nib is not yet properly defined 
at this stage. A procedure for actually making 
use of the previous formula is implemented in [5] . 
First of all, one matches the QCD propagator to 
its lattice HQET counterpart to obtain a relation 
involving the pole mass 

ml"'" = MB-£ + Sm + 0{l/mb). (2) 

*Talk presented by F. Di Ronzo 



This relation contains the quantity whose pcrtur- 
bative computation we are reporting on, i.e. the 
so called residual mass dm, a linearly divergent 
additive mass counterterm which is peculiar to a 
hard cut-off regularization scheme like the lattice. 
The pole mass can in turn be related to the MS 



mass rrib 



m 



MS 



/ \ pole 

mb[mb) = ml 



ri=0 



(3) 



the Dn coefficients with n < 2 being known from 
|5I7| . If we are now able to compute (ag is the 
lattice coupling) 



r!>0 

one can put everything together to get 

X„ 



(4) 



mb{mb) 



Mb ~ £ + Y.{a,r+^- 

oc 

A ^ TT 



(5) 



In this relation is a^ijfTb) (and the Xn get thus 
translated into X„). £ can be computed as the 
decay constant of the correlation function of two 
axial currents. Knowing 5m is crucial since it 
cancels both the divergence of £ and the renor- 
malon ambiguities which come from the matching 



between 



pole 



and 



,MS 



IH]. Since both cancela- 
tions take place in Perturbation Theory, one has 



1 



2 



to carefully assess to which accuracy one can con- 
trol the procedure. While Xq has been known 
for a long time, Xi was computed in while 



X 



(n/=0) 



was computed in and then also in 



[TT] . The knowledge of X. 



(»/=o) 



greatly improved 



the accuracy of the quenched determination of the 
6-quark mass. 

2. The computation 

jj-^"/-2) computed using the same strategy 
used for Xj"^"^^ ^U). We computed the Wilson 
loops W{R, T) for all the values of R and T up to 
16. This was done on a 32^ lattice using Numer- 
ical Stochastic Perturbation Theory. The com- 
putation was done for Wilson gauge and Wilson 
quarks action and the sea quarks masses were put 
to zero {i.e. the appropriate counterterms for the 
perturbative critical mass were plugged in). From 
Wilson loops and Creutz's ratios one can compute 

fW{R,T 



Vt{R)^ log 



1 



(6) 



V W{R,T) 
which in turn can yield the static potential via 
V{R)^YimVT{R). (7) 

1 — *oo 

The static potential is just the sum of the 
Coulomb potential (which can be described in 
terms of a potential coupling in whose definition 
the logarithmic divergencies^ are absorbed) and 
the (linearly divergent) mass counterterm we are 
interested in [T^ 

av(i?) 



V{R) ^25m - Cf- 



R 



(8) 



This means that the perturbative computation of 
the static potential can be read as the computa- 
tion of the matching between the lattice and the 
potential couplings, i.e. in the previous formula 



ay{R) — ao + ci{R) + C2{R) + 



(9) 



the coefficients ci (i?) and C2 (R) being dictated by 
A-parameters and /3-functions coefficients 

ci(i?) = 26ologi? + 26olog^ 
1^0 



^In Wilson loops also the so called corner logarithmic di- 
vergencies are present, but they cancel in the Creutz's 
ratios. 



and 



C2(i?) 



ci {Rf + 25i log R + 2bi log ^ 



r,(0) 



Since the matchings between both the potential 
and the MS couplings 13^ and the MS and the lat- 
tice couplings ri4] are known, the only unknown 
quantity in Eq. is Sm. Our procedure was to 



fit the coefficients X 



(n/=2) 



by matching our per- 



turbative computation of the static potential to 
Eq. |(SJ). An example of the fitting procedure can 
be seen in Fig. 1. The intervals in which fits were 
performed were such that: i? > 3, T > 12 and 
T > 2.5R {R is the mean value of R in the fit- 
ting interval; the fitting intervals themselves were 
from 3 up to 7 points long). We then choose in 
terms of ■ the errors we quote refer to the in- 
terval embraced by letting vary within a given 
interval. On top of that we could also inspect that 
the impact of lattice artifacts was under control 
(an handle is the value of known parameters en- 
tering the matching relations). 



3. Results 



We got Xo = 2.118(2), x[''f~^' = 10.56(4) 
(to be compared with the analytical values Xq = 



--2) 



2.118... and X[ 



("/=2) _ 



X. 



(»/=2) 



10.588...) and finally 

76.7(6). (10) 

Using the last result the authors of were able 
to repeat the analysis of their Uf = 2 unquenched 
data. The result for the 6-quark mass is 15 

rn^(mF)(™^) = (4.21 ± 0.03 ± 0.05 ± 0.04) GeV. 

At least three considerations are in order. First 
of all, by taking into account the value of Xj"^"^'* 
the central value moved from 4.26 to 4.21GeV and 
(what is more important) the last error, which 
is the one taking into account the indetermina- 
tions in the perturbative matching, got roughly 
halved. A second observation is that the new 
analysis yields different values for the quenched 
and unquenched results. Still, they are compati- 
ble within errors and in the end further investiga- 
tion is needed. The third point is the assessment 



3 




2 4 6 8 10 




^ ^ 



2 4 6 8 10 




2 4 6 8 10 

R 



Figure 1. The result of a fitting procedure as de- 
scribed in the text. Sohd Hues are the content of 
Eq. ^ once Sm has been fitted, while circles are 
the computed V^^\r), i.e. orders of Vt{R). 
In these figures T — 15. 



of sistematic errors involved in the procedure. In 
this approach one sticks to finite values of the 
lattice spacing. The finite lattice spacing depen- 
dence is not dramatic and it gets decreased by 
including the new term in the matching. On the 
other hand, the control on renormalon ambigui- 
ties seems firm. 

After the conference the analysis was refined and 
completed. The final version of the computation 
can be found in |16|. 
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